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ABSTRACT 

We use the Science Demonstration Phase data of the Herschel Virgo Cluster Survey to search for dust emission of early-type dwarf galaxies in the 
central regions of the Virgo Cluster as an alternative way of identifying the interstellar medium. We present the first possible far-infrared detection 
of cluster early-type dwarf galaxies: VCC781 and VCC951 are detected at the lOo" level in the SPIRE 250 yum image. Both detected galaxies have 
dust masses of the order of 10^ Mq and average dust temperatures » 20K. The detection rate (less than 1%) is quite high compared to the 1.7% 
detection rate for Hi emission, considering that dwarfs in the central regions are more Hi deficient. We conclude that the removal of interstellar 
dust from dwarf galaxies resulting from ram pressure stripping, harassment, or tidal effects must be as efficient as the removal of interstellar gas. 

Key words. Galaxies: dwarf; Galaxies: individual: VCC781; Galaxies: individual: VCC951; Galaxies: ISM; ISM: dust, extinction 
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1. Introduction 

Early-type dwarf galaxies (dEs) are the dominant morphological 
type in galaxy clusters. They were originally seen as a rather ho- 
mogeneous population of dwarf galaxies with an old stellar age, 
no features resembling any recent or ongoing star formation, 
and no indications of a significant interstellar medium (ISM). 
This viewpoint has changed radically in the past few years. Deep 
imaging observations of dEs have revealed a heterogeneous mor- 
phology. In particular, the population of dEs can be subdivided 
into nucleated and non-nucleated subclasses and several papers 
report evidence of disks, spiral structure, bars or star forma- 
tion (e.g., [Jerjen et al.||2000| IBarazza et al.|[20 0?, 'Geha et al~ 
2003] [&^m et al.||2003a|b( |De Rijcke et al . 2003 , Li sker et 
al.||2006a|b[ l. Kinematical studies demonstrate that dEs are not 



explosions) and external effects (ram-pressure stiipping, galaxy 
interactions, tidal effects) are thought to be able to expel the gas 



simple pressure-supported systems: some dEs seem to be rota- 
tionally supported ( Simien et al.|2002 Pedraz et al. 2002| Geha 
etal. |2003| |De Rijcke et al.|2003[|vaii Zee et a l. 2004; Toloba 
et al.|2009 1, whereas others show evidence of kinematically de- 
coupled cores (De Rijcke et al. 2004; Thomas e t al.|2006| l. 

Adding to this morphological and kinematical inhomogene- 
ity is the detection of a significant interstellar medium in vaii- 
ous dEs. Atomic and molecular gas was discovered for the first 
time in the Local Group dEs NGC 185 and NGC 205 using deep 
VLA a nd BIMA observadons d Young & Lo|1997| l. In the Virgo 
cluster, IConselice et al.| ( |2003t estimate a Hi detecti on rate of 
15% for dEs down to a 3cr Umit of 8 x 10^ Mq, while di Serego 
Ahghieri et al. ( 2007| l report a 1.7% Hi detection rate with a 3cr 
detection limit of 3.5 x 10^ Mq. Large amounts of cold inter- 
stellar matter are unexpected in dEs, as both internal (supernova 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



from the shallow potential on short timescales (Michielsen et al. 
i2004l Ivan Zee et a l. 2004, Roediger & Hensler 2005 , Boselh et 
al. 2008a b; Valck e'et al.]|'2008| ). That the Hi detected dEs are 
found preferentially near the edge of the cluster supports the 
idea that understanding environmental effects is crucial to con- 
straining the evolutionary history of dEs. Buyle et al.'(2005) and 
Bouchard et al. ( _2005) peiformed Hi observations in the out- 
skirts of the Fornax cluster and the Sculptor group, respectively, 
and confirmed the tendency of gas-deficient galaxies to be lo- 
cated near the center of the cluster. 

Continuum emission from interstellar dust is a promising al- 
ternative way to determine the ISM content of dEs. In paiticular, 
the ISM in the immediate suiToundings of M87 can be traced 
by dust emission, as strong radio continuum emission consid- 
erably reduces the ALFALFA Hi detection sensitivity within 1° 
of M87 (Giovanelh et al.|2007"l ). To date, the Andromeda satel- 
lites NGC 205 and NGC 185 are the only dEs that have been de - 
tected in the far-infrared ( |Haas|1998|rMarIeau et al.]2006l|2009] l. 
In paiticular, no dust emission has been detected from cluster 
dEs. While extinction features in optical images are indicative 
of dust in at least some cluster dEs ( |Ferrarese et al. 2006 ; Lisker| 
|et al.||2006a|b| l, the direct detection of dust emission has been 
hampered by the limited resolution and sensitivity of the pre- 
vious generation of far-infrared instrumentation. The advent of 
the Herschel Space Observatory ( Pilbratt et al.|2010| l offers new 
possibilities for mapping the ISM in dEs. In this paper, we report 
on our search for far-infrared emission from dEs in the Virgo 
Cluster, based on the Science Demonstration Phase (SDP) ob- 
servations of the Herschel Virgo Cluster Survey (HeViCSfM). In 
Sect. [2] we describe the observations, data reduction, and sample 



More details on HeViCS can be found on http://www.hevics.org 
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selection. In Sect. [3] we present our detections and in Sect. |4]we 
discuss the results of a stacking analysis of the non-detections. 
Sect.lSlgives our summary. 

2. Observations and sample selection 

The central 4x4 deg^ region of the Virgo Cluster was ob- 
served by Herschel on 29 November 2009 as part of the HeViCS 
SDP observations. The observations were performed in paral- 
lel scan map mode, which means that both PACS ( [Poghtsch et 
[ar] [20T0| and SPIRE ( Griffin et al.pOlO. ) data were coUected. 
The scan speed was 60"/s and the field was covered in nomi- 
nal and orthogonal directions to minimize 1 // noise. The PACS 
and SPIRE data were reduced with HIPE, with reduction scripts 
based on the standard reduction pipelines. For more details of the 
HeViCS SDP data reduction, we refer to Davies et al., ( ,201 0). We 
used the HeViCS SDP observations to search for dust emission 
from dEs in the central regions of the Virgo Cluster. Selecting the 
galaxies with morphological type dE and dSO in the Goldmine 
catalogue ( |Gavazzi et al. |20()3] l resulted in an initial sample of 
370 dEs. After the removal of 16 galaxies classified as back- 
ground sources and the exclusion of 115 sources with an addi- 
tional optical source within a 6" radius, our final sample con- 
sisted of 239 dEs. Blind aperture photometry was applied to the 
PACS and SPIRE maps based on the positions of each of these 
239 galaxies using Source Extractor (jBertin & Arnouts|1996i. 



3. Detected dEs 

From this sample of 239 sources, 1 1 dEs were detected above 
3cr in the SPIRE 250 yum image. Nevertheless, we only report 
the clear detection of 2 dEs,VCC781 and VCC951, exclud- 
ing all sources raising any kind of doubt. As such, VCC 1502 
and VCC 788 were ruled out because their FIR location dif- 
fers more than FWHM/2 from their optical position. VCC 832 
was excluded because SDSS images detect a background source 
with substructure in the residual image, possibly indicating 
the presence of dust. Four other objects (VCC 752, VCC 815, 
VCC 1272, VCC 1512) fall within a crowded field of sources, 
enhancing the chances of their being background sources. For 
instance, Binggeli et al.| (11985) reports the detection of a blue 
compact background galaxy, only 4" to the west of the nucleus 
of VCC 815. Furthermore, VCC 752 and VCC 1272 appear to 
have SEDs that do not correspond to a single grey body. This 
failure in fitting a SED might either be caused by poor flux de- 
termination (the error bars are large compared to the low number 
of counts in the detection) or be a true indication that both detec- 
tions originate in sources other than dEs. The other two sources 
VCC 1578 and VCC 767 are clearly detected in the SPIRE 250 
jum image at 4.8 and Act, respectively, but their high FIR flux 
relative to their faint optical appearance (both objects are low 
surface brightness [LSB] galaxies with B-band magnitudes of 
the order ~2I mag) aroused suspicion. 

Dust has already been detected in LSBs ( Hinz et aL]|2007 i, 
but mainly in those with large amounts of star formation respon- 
sible for dust heating. A lack of evidence of recent star forma- 
tion in VCC 767 and VCC 1578 strongly argues against the as- 
sociation of the 250 //m emission with the corresponding dEs. 
Furthermore, we can identify a background source in the imme- 
diate surroundings of VCC 1578, which possibly contaminates 
the SPIRE 250 //m emission. 

Nonetheless, these questionable detections may still be 
caused by dust emission from dEs. Deeper fields in the future 



might reveal dust emission in PACS images, where the higher 
spatial resolution may be able to rule out possible contamination 
by background sources. 

On the other hand, VCC 781 and VCC 951 are detected un- 
ambiguously in the SPIRE 250 yum image at 9.8 and ll.lcr, re- 
spectively (see Fig. fT). Adopting a fair 2cr detection cutoff in 
the other filters, VCC 781 is also detected in the PACS 100 yifm, 
PACS 160 yum, SPIRE 350 yum, and SPIRE 500 jum bands. 
VCC 951 has additional detections only in both SPIRE 350 //m 
and 500 jiva images. Unfortunately, it is impossible to measure 
fluxes for VCC 951 at 350 jum and 500 jum because VCC 951 
and its neighbouring source, SDSS J122655.56-t-l 14000.5, are 
not resolved at these longer wavelengths. 

Since the SPIRE 250 yum image contains a vast number of 
sources, we need to address the possibility of false detections 
due to background sources. To perform this check accurately, we 
search all relevant images/catalogs (SDSS, Goldmine, 2MASS, 
UKIDSS, Spitzer) for sources within 18.1 "(FWHM of SPIRE 
250 yum) of the dwarf positions. For VCC 78 1 , this identifies four 
SDSS sources within this PSF region that may be contaminating 
the dwarf emission. The superior resolution in the PACS 160 yum 
image discards three of these optically identified objects, while 
the remaining source is rejected as a possible background candi- 
date for the IR emission, based on the 24 /im image (data from 
Fadda et al. in prep.), which has an even higher spatial resolu- 
tion, and is clearly extended. Hence, the coincidence of the po- 
sition of detections at various wavelengths makes us confident 
that the IR emission originates in VCC 781. 

The SDSS DR7 catalog provides us with 5 detections within 
the PSF region of VCC 951, three of which are very likely to be 
artifacts in the source extraction, since we detected them in nei- 
ther the SDSS nor the UKIDSS images during our visual inspec- 
tion. The remaining two, morphologically classified as stellar- 
like objects, are located at distances of 8" and 10", and are also 
detected in the UKIDSS images. The 24 yum image (see Fig. [Til 
shows some evidence of emission from both VCC 951 and one 
of the two aforementioned sources, which is spatially resolved 
at these wavelengths. Based on these data, we cannot definitely 
conclude that the IR emission detected by SPIRE comes from 
our dwarf galaxy, but deeper optical and NIR observations may 
help us to address the issue. 

We conclude that, while our analysis is unable to definitely 
exclude a possible contamination from background sources, the 
weight of evidence is in favour of a true detection of dust in 
emission from dEs, and in particular from VCC 781. 



These detected dEs have remarkable morphologies. Binggeli 
[eTaL] ( fT985] l classifies VCC781 as dS03(5),N and VCC951 
as dS0(2),N or dE2 pec,N. These galaxies are also found by 
[Lisker et al.| (|2006all to harbor central substructures other than 
a disk, while [Lisker et a l. (2006b) identified blue central colors 
in VCC 781 and VCC 951. Figure [TJpresents the g - i color im- 
ages in which VCC 781 and VCC 951, respectively, exhibit an 
obvious gradient in their radial g - i color profiles, strengthen- 
ing the assumption of recent star formation in the central regions 
(see also |Lisker et a l. 2006b). 

According to the classification criterion (3 < FUV-H < 6) 
adopted in [Bosem et aL] (p008a), VCC 781 and VCC 951 can 
also be classified as possible transition objects. Although the 
FUV-H colours of 6.6 and 6.9 mag for VCC 781 and VCC 951, 
respectively, do not satisfy this relation, both galaxies are clearly 
located at the blue end of the dE galaxies in Boselli et al. 
(2008a I. That SDSS nuclear spectra of both galaxies exhibit 
deep Hydrogen absorption lines (EW[H5] > 4 A) indicates that 
they are in a post-starburst phase. This, and their 24 yum emis- 
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Fig. 1. All images are on the same scale. The FWHM of the beam is denoted in the MIPS 24 //m and SPIRE 250 yum images, g-band 
contours are overlaid on the SPIRE 250 ;um image. The dwarfs position is indicated as a star, SDSS background sources as open 
circles, and the neighbouring source of VCC951 as a cross. 



sion, which is indeed centered on the optical nuclei and concen- 
trated (~ 6" for VCC781 and point-like for VCC951), argues 
in favour of a connection between dust emission and a recent 
episode of star-formation. 

To estimate the dust masses, we determined fluxes in bands 
where we had detections > 2cr. Initial apertures were fixed in 
the SPIRE 250 //m images based on the flattening of the growth 
curve. These apertures were subsequently adjusted to the pixel 
scale in other bands such that all apertures cover the same phys- 
ical area. We first determined a representative dust temperature 
using the assumption that the dust is in thermal equilibrium with 
the interstellar radiation field. We used the Monte Carlo code 
SKIRT, which was initially developed to investigate the effects 
of dust extinction on the photometry and kinematics of galax- 
ies (Baes et al.||2003| l, but evolved into a flexible tool that can 
model the absorption, scattering, and thermal emission of cir- 



cumstellar discs and dusty galaxies (e.g Vidal & Baes 2007 



Baes et al. 2010 1. The stellar body of each galaxy was repre- 



sented as an exponential profile, with parameters taken from the 
Goldmine database. The dust was assumed to have the same dis- 
tribution as the stars. For the intrinsic SED of the model, we 
used the elliptical galaxy template SED from the PEGASE li- 
brary (Fioc & Rocca-Volmerange| |1997 1. Adding these ingredi- 
ents, we find a representative dust equilibrium temperature of 
20.7 K and 19.4 K for VCC781 and VCC951, respectively. 
Relying on these temperatures, we compute the dust masses ac- 
cording to 

Md=^^5^, (1) 

''■250 fi250(7'd) 



where D - 16.5 Mpc is the distance to the Virgo Cluster (Mei 
et al. 2007 i, K2so - 0.4 m^ kg"' is the dust absorption coefficient 
( |Draine & Li|200l"] l, and B25o(rd) is the Planck function for the 
modeled temperature T^. We obtain dust masses of 1.85x10^ Mq 
and 1.28 x 10^^ Mq for VCC781 and VCC951, respectively. 

As a check for consistency with datapoints other than the 
SPIRE 250 yum flux, we determine a second estimate of the tem- 
perature and dust mass based on grey body fitting. A single grey 
body fit to the 5 PACS and SPIRE fluxes gives Tj = 19.5 K and 
Md = 2.24 X 10^ Mq for VCC 781 (see Fig.|2]i. These results are 
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Fig. 2. The best-fitting grey body for VCC 78 1 with Tj = 19.5 K 

and dust mass Md = 2.24 x 10^ Mq. 



fully consistent with the temperature and dust mass of VCC 781 
derived above. The MIPS 70 and 160 jim datapoints also satisfy 
the fitted greybody curve well (see Fig.l2|. 

Neither of the two dwarf galaxies is detected in Hi. An Hi 
mass upp er limit was set to 2.3 x 10^ M© for VCC 78 1 (|Gavazziet| 



al. 2003ll and 8.0 x lO*" M^ for VCC 95 1 dConselice et al. 2003 



Combining the estimated dust masses with these Hi upper limits, 
and neglecting molecular gas, we find Hi gas-to-dust upper limits 
of 124.3 and 62.5, respectively. These values are much lower 
than the expected canonical gas-to-dust ratio, possibly indicating 
that there is a dependence on other factors such as the metallicity 
and the environment. 



4. Stacking analysis of non-detections 

Apart from the two detections described in the previous section, 
we performed a stacking analysis of non-detected sources. In a 
first step we masked all bright sources in the surroundings of the 
dwarf galaxy (we prevented the inner pixel from being masked). 
We enlarged these masks using a boxcar with smoothing ker- 
nel of 3 pixels to cover the whole source. The final stacking al- 
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gorithm weighs the images with the inverse of the background 
variance to acquire the mean value in each pixel (of all unmasked 
pixels). 

Stacking the SPIRE 250 //m images of all 227 non-detected 
early-type dwarf galaxies results in no detection. When taking 
apertures of size equal to the SPIRE 250 jum beam at random 
positions and calculating the standard deviation in their distri- 
bution, we derive an estimated background rms of 0.38 mJy. 
Assuming an average dust temperature of Tj - 20K, we con- 
vert this noise estimate to a 3cr dust mass detection limit of 
Md = 2.44 X 10^ Mq. 

Since the noise measured locally has an average value of rms 
= 5 mJy/beam, the stacking is not significantly affected by con- 
fusion noise (the rms almost scales as Gaussian noise, 5/ V227 
= 0.33 * 0.38). The complete survey of 4 cross-scans will reach 
a Icr sensitivity of 4.72 mJy on small spatial scales, considering 
both instrumental and confusion noise (the latter is determined to 
be 4 mJy at SPIRE 250 //m). Assuming that the stacking proce- 
dure is only marginally affected by confusion noise, this value 
translates into a 3cr detection limit of Md = 2.0 x lO"* M©. 
Comparing both values, we conclude that the stacked image 
nearly reaches the sensitivity limit, as the detection limit will 
only improve by a factor of 1 .2 after the complete survey. 

5. Discussion and conclusion 

We have reported on our search for far-infrared dust emission in 
a sample of 239 dEs at the centre of the Virgo Cluster With a 
detection limit of about lO'' M©, we have detected two objects at 
lOcr that have a high probability of being true detections of dEs. 
These two dwarf galaxies, VCC 781 and VCC951, are the first 
detections of dust emission in dEs, apart from the Andromeda 
satellites NGC205 and NGC 185. Applying a radiative transfer 
model to determine the equilibrium dust temperature in their in- 
terstellar radiation field, we determined temperatures of the or- 
der of 20 K in both galaxies. We estimate dust masses of the 
order Md = 10^ M© in both dEs. 

In addition to independent detections of both atomic and 
molecular gas in dEs, the detection of dust emission here con- 
firms the hypothesis that at least some dEs are transition objects 
gradually evolving from late-type to early-type, while falling 
into the cluster 

The detection of an interstellar dust medium in two dEs is 
rather unexpected, given the efficiency with which the ISM is 
removed from dwarf galaxies. Boselli et al.| ([2(308a) found that 
gas removal by ram pressure stripping is extremely efficient in 
dwarfs galaxies, there being a decrease by ~2 orders of mag- 
nitude in ~150 Myr of time. In agreement with this rapid mor- 
phological transition, only a small number (10 - 16%) of dwarf 
galaxies exhibiting properties intermediate between star-forming 
and spheroidal galaxies are present in the Virgo Cluster. 

di Serego Alighieri et al. ( 2007| l find an Hi detection rate 
in the Virgo cluster of only 1.7% to a 3cr mass limit Mhi - 
3.5 X 10^ Mq. Adopting a canonical gas-to-dust ratio of 1000 
(which might still be a conservative lower limit as dEs in Virgo 
are expected to have metallicities significantly lower than ob- 
served in our own Milky Way), we recover a 3cr dust mass limit 
of S.SxIO'^Mq, which is more or less the 3cr dust mass detection 
limit in the HeViCS SDP field. Our detection rate is also compat- 
ible with the 15% found by Conselice et al. ( |2003[ l. Their Hi up- 
per limit corresponds to a 3cr dust mass limit of 8 X 10"^ Mq, well 
below the SDP detection limit. When also taking into account 
the position in the cluster, we consider our detection rate of less 



than 1 % surprisingly high. Indeed, Hi detected dEs are prefer- 
entially located in the outskirts of the cluster and may not have 
been dramatically affected so far by the intracluster medium and 
interactio ns with other galaxies ([Consehce et al . 2003 ; B uyle et| 
[aLl|2005| l. |di Serego Alighieri et al.| ( |2007| l also report that the 
majority of the detected Hi-sources are located near the edges of 
the Virgo cluster If gas and dust are tightly coupled and dust is 
as easily removed from the galaxy as gas by ram pressure strip- 
ping, harassment, or tidal effects, we expect that the dEs in the 
central regions of the Virgo Cluster will also be the most dust 
deficient. Clear evidence of efficient dust stripping in the Virgo 
Cluster is seen for the first time in Hi-deficient spiral galaxies 
(Cortese et al..2010) . Future HeViCS observations will be able 
to provide an answer to this question - the HeViCS survey will 
not only go deeper (each field will be covered 4 times), but also 
wider (the total area covered will be ~64 deg^). A comparison 
of the far-infrared detection rate with the position in the cluster 
will tell us whether environmental effects have any influence on 
the detection rate. 
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